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ABSTRACT 


A  SIMULINK  towing  model  is  developed  from  the  surge,  sway,  and  yaw 
equations  of  motion  in  order  to  study  the  horizontal  maneuverability  of  vessels  in  a  semi¬ 
rigid  towing  operation.  This  analysis  is  eonducted  in  order  to  validate  rigid-eonneetion 
towing  and  to  give  insight  into  the  design  of  the  tow  conneetor.  The  connection  is 
modeled  as  a  linear  spring  and  the  maneuverability  of  the  vessels  is  studied  as  the 
stiffness  is  varied  from  conditions  of  semi  to  completely  rigid.  This  study  is  based  on 
two  Swath  hull  vessels,  the  SLICE  and  KAIMALINO,  towing  in  close  proximity  under 
calm  water  conditions. 
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I.  INTRODUCTION 


As  part  of  a  continuing  program  to  advance  in  the  swath  hull  teehnology,  this 
paper  looks  into  the  feasibility  of  towing  operation  between  the  Slice  and  Kaimalino. 
The  Sliee  is  an  advaneed  high-speed  vessel  that  is  emerging  as  a  ship  with  a  very  wide 
range  of  open  oeean  eapabilities.  The  Offiee  of  Naval  Researeh  desires  to  further 
showease  the  advantage  of  multi-hull  design  by  demonstrating  the  maneuverability  of  the 
Slice  and  Kaimalino  during  towing  operations.  Although  numerous  investigations  have 
been  condueted  on  optimizing  towing  connection  points,  towline  length,  and  tension  of 
towline,  this  study  will  be  eoneerned  with  elose  proximity  towing  under  a  near  rigid  tow 
eonneetor.  Furthermore,  most  previous  studies  model  the  dynamics  of  the  trailing  ship 
only  while  the  leading  ship  is  treated  as  a  point  mass.  This  thesis  models  the  eoupled 
maneuvering  motions  of  the  two  vessels. 


A.  KAIMALINO 

The  SSP  Kaimalino  was  the  world’s  first  high  performanee  open  oeean  Swath 
ship.  It  eonsists  of  two  parallel  torpedo-like  hulls.  Attaehed  to  the  hulls  are  two 
streamlined  struts.  The  struts  extend  above  the  water  surfaee  and  support  the  main  body. 
The  Kaimalino  also  has  stabilizing  fins  attaehed  near  the  aft  end  of  eaeh  hull.  Figure  1 
shows  a  profile  of  the  SSP  Kaimalino.  The  performanee  features  of  the  Kaimalino  are 
greatly  reduced  motions  with  sustained  speed  even  in  high  sea  states,  lower 
hydrodynamie  drag  and  redueed  power  requirements  at  moderate  to  high  speeds,  and 
excellent  eourse  keeping  eharaeteristies  at  all  sea  headings. 
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Figure  1.  SSP  KAIMALINO  General  Arrangements  from  Referenee  [8]. 
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Length 

88.6  ft 

Beam 

45.9  ft 

Height 

31.8  ft 

Displacement 

217  tons 

Payload 

50  tons 

Max  Speed 

25  kt 

Operational  Speed 

13  kt 

Range 

400  n  miles 

Table  1.  KAIMALINO  characteristics. 


B,  SLICE 

The  Slice  is  a  high  speed  variant  of  the  Swath  technology.  It  has  4  underwater 
hulls  instead  of  two.  Attached  to  each  hull  is  a  strut  that  extends  up  to  support  the  main 
body.  Figure  2  and  3  show  the  profde  view  and  the  stem  view  of  the  Slice  respectively. 
The  Slice  produces  higher  speeds  for  constant  horsepower  when  compared  to 
conventional  Swath  hulls  and  still  maintains  the  good  sea-keeping  of  Swath  technology. 
The  key  innovation  is  reduction  of  wavemaking  drag  accomplished  by  using  short  stmts, 
submerged  hulls,  and  speeds  well  beyond  the  critical  “hump”  on  the  Froude  resistance 
curve  [1]. 
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Figure  2.  Profile  view  of  SLICE  vessel  from  Referenee  [10]. 


Figure  3.  Stem  view  of  SLICE  vessel  from  Referenee  [10]. 
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Table  2.  SLICE  characteristics. 
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II.  MANUEVERING  ANALYSIS 


A,  EQUATIONS  OF  MOTION  IN  THE  HORIZONTAL  PLANE 

For  the  case  of  horizontal  plane  maneuvering,  the  vertical  plane  motions  are 
ignored,  leaving  only  the  surge,  sway  and  yaw  equations  of  motions  (EOM)  for  the 
equation  set.  These  equations  will  be  developed  with  the  assumptions  that  xg  and  yo  are 
zero,  and  vehicle  symmetry  for  its  inertial  properties.  Also  in  the  derivation  of  the  EOM, 
two  sets  of  axes  are  used,  the  inertial  frame  and  the  body  fixed  frame.  The  inertial  frame 
is  relative  to  the  earth,  ignoring  the  earth’s  rotational  rate  versus  the  vessel’s  rotational 
rate.  Newton’s  Eaw  is  applied  in  the  inertial  reference  frame.  The  body  fixed  frame  is 
relative  to  the  ship’s  center  of  gravity  having  the  positive  x-axis  forward  of  the  ship, 
starboard  for  positive  y-axis,  and  down  for  positive  z-axis.  Applying  Newton’s  law  in  the 
inertial  reference  frame  results  in  the  following  equations  for  surge,  sway,  and  yaw. 


(1) 

«To  =  Yo 

(2) 

II 

(3) 

where  Xq  and  Yq  are  the  total  forces  in  respective  directions  and  N  is  the  turning 
moment  about  the  z-axis.  Detailed  in  Reference  [4]  are  the  components  that  usually 
make  up  the  total  forces  and  moments.  These  total  forces  are  a  combination  of  fluid 
forces  acting  on  the  main  body,  control  surface  forces  (i.e.  rudders,  fins,  bow  planes), 
environmental  forces,  and  propulsion  forces.  Eooking  at  the  diagram  below,  it  is  easily 
shown  that  a  relationship  can  be  developed  for  transformation  from  inertial  reference 
frame  to  the  body  fixed  frame. 
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Figure  4.  Relationship  of  position  between  two  reference  frames. 


Therefore  we  transform  coordinate  systems  using 


X  =  Xq  cos  y/+Y^  sin  y/  (4) 

Y  =  cos  y/  -  sin  y/  (5) 

and  for  similar  velocity  diagrams 

x  =  ucosy/ -vsiny/  (6) 

y  =  usmy/  +  vcosy/  (7) 


The  rate  of  yaw  is  given  by  y/  =r.  Then  differentiating  the  above  equations  with 
respect  to  time  results  in 

x  =  ucosy/ -uy/siny/ -vsiny/ -vy/cosy/  (8) 

y  =  ucosy/ -uy/siny/ -vsiny/ -vy/cosy/  (9) 
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Substituting  these  results  into  Equation  (1)  through  (3)  then  into  Equations  (4)  and 
(5)  results  in  the  following, 


mil  —  mvr  =  X 


mv  +  mur  =  Y 
I,r  =  N 


(10) 

(11) 

(12) 


where  the  u  is  the  surge  veloeity,  v  is  the  sway  velocity  and  r  is  the  yaw  velocity. 


Now  expressions  for  the  right  hand  side  (the  forces)  need  to  be  developed.  The 
environmental  forces  will  be  neglected  for  this  study  and  the  propulsion  force  will  be 
designated  only  as  Xprop.  Neglecting  any  side  propulsion  forces  resulting  from  torque,  the 
total  propulsion  force  acts  only  in  the  longitudinal  direction  for  the  surface  vessel.  Eor 
the  fluid  and  rudder  forces,  the  nonlinear  components  are  complicated  and  expressions  to 
approximate  those  components  are  still  in  development.  The  fluid  forces  will  be 
linearized  using  a  Taylor  series  expansion  of  forces  as  steady  state  functions  in  order  to 
produce  body  forces  represented  as  hydrodynamic  coefficients.  The  hydrodynamic 
forces  are  functions  of  vessel  velocities  and  accelerations.  Expressions  will  be 
developed  for  small  deviations  from  a  nominal  condition  of  straight  line  transit  at  a 
constant  speed  El. 


Yp  ={YF)o  +  ^{u-UQ)  +  —{v-v^)  +  —{r-rQ)  +  —  {u-UQ)  + 

au  av  or  au 

ov  or 


(13) 


where  the  subscript  E  is  used  for  fluid  force. 


Eor  the  nominal  condition  stated  above,  the  initial  conditions  are  as  follows. 
Uq  =U 

Vo  =  Wo  =  Vq  =  To  =  fg  =  0 
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Since  the  ship  is  aeting  in  straight  line  motion  there  is  no  initial  transverse  foree. 
Also  beeause  of  symmetry,  changes  in  the  longitudinal  veloeity  and  aeeeleration  produce 
no  transverse  foree,  dYjdu  anddYjdu  are  zero. 


37.,,  37.,,  37.,.,  37.,., 
Y.  =  ^(v)  +  ^(r)  +  ^(v)  +  ^(r) 


3v 


3r 


3v 


3f 


(14) 


Similar  logie  leads  to  the  resulting  expressions  for  surge  and  yaw. 


N, 


du 


(u  —  1/)  + 


du 


dN,  dN,  .  dN,  dN, 


dv 


-v  +  - 


dv 


-v  +  - 


dr 


-r  +  - 


dr 


(15) 


(16) 


The  following  notation  for  the  derivatives  will  be  used  for  the  remainder  of  this 
paper;  3X/3u  is  Xu,  3Y/3u  is  Yu,  dX/dii  is  Xi,  ete...  These  are  the  hydrodynamie 
derivatives. 

Assuming  that  the  rudder  forees  and  moments  on  the  ship  are  funetions  of  delta 
only,  they  too  can  be  linearized.  The  rudder  deflection  angle  is  represented  by  5.  A 
positive  5  means  vessel  is  turning  to  port  and  negative  5  is  for  turn  to  starboard. 
Expressions  for  the  rudder  forees  and  moments  are  presented  below. 


X,=X,(cJ  =  0)  +  ^(^  =  0 
do 

(17) 

(18) 

N,=M,(S  =  0)  +  ^S  =  N,S 

(19) 

The  subseript  R  is  used  for  rudder  forees. 

Substituting  back  into  equations  (10)  through  (12),  we  get  the  following  EOMs 
for  surge,  sway  and  yaw. 
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im-X„)i  =  X,(u-U)  +  +  X,S 

(m  -  Y.)v  -  Y-r  =  Y^v  +  (7^  -  mU)r  +  YgS 


(20) 


(21) 

(4  -N,)r-  N,v  =  N^v  +  N/  +  NgS  (22) 

B.  EQUATIONS  OF  MOTION  FOR  THE  SLICE  AND  KAIMALINO 

In  developing  the  EOMs  for  the  towing  model,  another  force  had  to  be  added. 
This  force  results  from  the  tension  in  the  towing  connection.  Using  only  the  forces  in  the 
horizontal  plane,  this  force  is  expressed  as  some  variable  spring  constant  K  times  the 
difference  between  the  SLICE  and  Kaimalino  positions.  The  towing  operation  to  be 
modeled  is  shown  in  Eigure  5. 


Eigure  5.  Diagram  of  tension  forces  resulting  from  vessels  under  tow. 
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T,=kMc,x-^c,2) 


(23) 


Ty  =  K{y,,-y,^,)  (24) 

subscript  1  represents  the  SLICE  and  subseript  2  is  for  the  Kaimalino.  This  linear 
relationship  expressing  the  tension  forees  eouples  the  SLICE  and  Kaimalino  EOMs.  Ky 
and  Kx  are  taking  at  various  values  in  order  to  analyze  the  maneuverability  of  the  vessels 
in  towing  operation.  Below  are  the  eomplete  EOMs  for  the  SLICE  and  Kaimalino. 


{m\-  Xui)u\-m\riv\  =  -R\  +  Xprop-TxCosy/\  +  Tysmy/\  (25) 

(mi  -  Tvi)vi  -  Kifi  +  miriMi  =  Yvxvx  +  Ynrx-TyCosif/x-Txsmy/x  +  Ysd  (26) 

(/zi  -  Nn)r\ -  Nmvx  =  Nvm  +  Nnr\  +  Xp\{TxS,my/\  +  Tycosy/i)  +  Nsd  (27) 

{mi  -  Xu2)u2  -  mirivi  =  -Ri  +  TxCOSif/i  +  Tysmy/i  (28) 

{m2-Yv2)v2-Yr2r2  + miriui  =  Yvivi  +  Yr2r2  +  Tycosi/^i  -  Txsini/^i  (29) 

{Iz2- Nr2)r2- NnVl  =  Nv2V2  +  Nr2r2  +  Xp2{TyCOSy/2-TxSmy/2)  (30) 


where  y/\  and  y/i  are  the  yaw  angles  of  the  SLICE  and  Kaimalino  respeetively. 
Ri  and  R2,  the  resistanee  of  the  vessels  moving  through  body  of  water,  are  funetions  of 
the  ships  speeds  through  the  water.  Eor  this  analysis,  it  will  be  assumed  that  the 
propulsion  from  the  SLICE  adjusts  in  order  to  maintain  a  eonstant  speed  for  U. 
Therefore,  the  surge  equation  for  the  SLICE  will  be  ignored  and  ui  ean  be  taken  as  a 
eonstant. 

The  following  equations  are  formed  from  transformations  between  eoordinate 
systems  and  provide  the  neeessary  relationships  in  order  to  solve  the  total  systems  of 
equations. 


Xj  =  Mj  eos  y/^  -  Vj  sin  y/^ 


(31a) 
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}^  =  Mj  sin  +  Vj  cos  (31b) 

Xj  =  U2  cos  y/2  -  V2  sin  XI/2  (32a) 

Y2  =  U2  sin  y/2  +  Vj  cos  y/2  (32b) 

also 

Wx  =  (33a) 

y/2  =  ^2  (33b) 


C.  SIMULINK  TOWING  MODEL 

The  towing  model  is  built  using  the  Matlab  Simulink  toolbox.  The  dynamie 
responses  will  also  be  analyzed  using  Simulink.  Simulink  allows  for  eontinous  time 
modeling  of  both  linear  and  nonlinear  models.  Eaeh  of  the  EOMs  developed  above  are 
modeled  in  Simulink  as  superbloeks.  Each  of  these  equations  will  be  eoupled  together  by 
the  expressions  that  were  built  for  the  tension  in  the  towing  connector.  One  additional 
transformation  needs  to  be  developed  sinee  we  are  interesting  in  analyzing  the 
maneuvering  of  the  vessels  due  to  towing.  The  equations  of  motions  are  developed  for 
the  X  and  y  values  being  at  the  eenter  of  gravity  (assumed  midships).  The  equations 
below  are  developed  in  order  to  reference  ship  motions  to  the  point  of  tow  conneetion 
rather  than  eenter  of  gravity. 


Xc,x  =  X,-^cosy/, 


V  ■ 

Tc.i 


X,2  =^2+yCos^^2 


Tc,2  =  1^2+^8111  ?^2 


(34a) 

(34b) 

(34e) 

(34d) 
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The  subscript  c  is  used  for  the  tow  connection  position.  Both  the  tension 
expression  and  above  expressions  for  tow  connector  positions  are  also  modeled  with 
superblocks. 

The  Simulink  towing  model  is  shown  in  Figure  6.  The  diagram  shows  how  the 
Slice  sway  and  yaw  EOMs  are  coupled  together  by  the  state  variables  vi  and  ri.  The 
Kaimalino  surge,  sway  and  yaw  EOMs  are  coupled  by  the  state  variables  U2,  V2,  and  r2. 
Each  EOM  gives  an  output  state  that  feeds  into  coordinate  transformation  equation  and 
results  in  values  of  X  and  Y.  The  linear  towing  model  couples  the  motions  of  the  Slice 
with  that  of  the  Kaimalino.  Simulink  integrator  blocks  are  used  to  get  the  state  of  certain 
variables.  Internal  initial  conditions  are  used  for  the  integrator  blocks  that  allow  for 
setting  starting  values  such  as  position  and  orientation. 

TOWINQ  MODEL  FOR  SUCE  AND  KAiMAUNO 


Eigure  6.  Diagram  of  Simulink  Towing  Model  for  Slice  and  Kaimalino. 
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The  sway  equation  of  motion  for  the  Slice  and  Kaimalino  is  shown  in  Figures  7 
and  8.  The  diagrams  are  similar  with  the  exception  of  rudder  force  term  that  is  only  in 
the  Slice  diagram.  The  longitudinal  velocity  (u),  yaw  velocity  (r),  and  yaw  acceleration 
(f  )  are  input  state  variables.  The  lateral  velocity  is  an  output  state  variable  and  is  fed 
back  into  the  equation  using  an  internal  loop.  Each  one  of  the  state  variables  is 
multiplied  by  respective  force  or  hydrodynamic  derivative  using  gain  blocks.  The  gain 
blocks  also  include  the  dividing  through  of  m-Y-  term.  The  variable  used  for  gain 

blocks  are  defined  in  an  m-file  that  initializes  the  Simulink  model.  The  tension  terms  are 
like  forcing  functions  and  are  added  to  the  system.  Each  one  of  the  tension  terms  is 
multiplied  by  the  gain  m-Y..  The  yaw  angle  {y/),  tension  (Ty,  Tx),  u  and  r  are  provided 

as  inputs  from  other  superblocks.  The  Slice  uses  a  constant  u,  whereas  the  Kaimalino 
uses  the  input  from  the  surge  equation.  The  equation  and  diagram  are  designed  to 
produce  the  single  output  of  the  lateral  acceleration  v . 


SUCE  SWAY  EQUATION  OF  MOTION 


Eigure  7.  Diagram  of  Superblock  for  Slice  Sway  EOM. 

15 


THE  KAIMALINO  SWAY  EQUATION  OF  MOTION 


Txslnpsl2 


Figure  8.  Diagram  of  Superblock  for  Kaimalino  Sway  EOM. 


The  yaw  equations  of  motion  for  the  Slice  and  Kaimalino  are  shown  in  Figures  9 
and  10.  The  process  for  building  the  superblock  is  similar  to  that  used  for  the  sway 
equation  of  motion.  The  use  of  gain  blocks  are  used  again  here  to  represent  the 
hydrodynamic  coefficients  divided  through  by  -  N-  term.  The  single  output  for  the 
yaw  superblocks  is  yaw  acceleration  r  and  the  yaw  velocity  r  is  fed  back  into  the 
subsystem  using  an  internal  loop. 
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THE  SLICE  YAW  EQUATION  OF  MOTION 


Figure  9.  Diagram  of  Superblock  for  Slice  Yaw  EOM. 

THE  KAIMALINO  YAW  EQUATION  OF  MOTION 


Figure  10.  Diagram  of  Superblock  for  Kaimalino  Yaw  EOM. 
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Figure  1 1  shows  the  superbloek  that  represents  the  Kaimalino  surge  equation  of 
motion.  For  this  superbloek  a  look-up  table  was  used  to  model  the  resistanee  term  in  the 
surge  equation.  The  lookup  table  matehes  an  input  to  an  output  using  linear  interpolation 
built  into  the  eode.  The  input  for  this  resistanee  is  the  longitudinal  speed  of  the 
Kaimalino.  A  graph  representing  the  input  and  output  data  is  plotted  on  the  lookup  table 
ieon.  Simulink  lookup  table  ean  be  used  for  n-dimensions,  therefore  other  resistanee 
forees  ean  be  modeled  sueh  as  wind  and  eurrent.  The  other  terms  and  the  surge  equation 
are  modeled  similar  to  the  gain  bloeks  used  in  the  sway  and  yaw  equations  of  motion. 
The  single  output  for  this  superbloek  is  the  Kaimalino ’s  longitudinal  veloeity  u. 


THE  KAIMALINO  YAW  EQUATION  OF  MOTION 


Figure  1 1 .  Diagram  of  Superbloek  for  Kaimalino  Surge  EOM. 


All  the  variable  inputs  from  the  equations  of  motions  for  the  Sliee  and  Kaimalino 
are  used  in  the  ealculation  of  X  and  Y.  Figure  12  shows  this  superbloek.  The 
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transformation  Equations  (31a)  through  (31b)  are  modeled  here.  The  initial  eonditions  to 
set  up  the  vessel  spaeing  and  length  of  tow  are  entered  using  the  integrator  bloeks  for  the 
X  and  Y  states.  The  values  are  in  terms  of  ship  lengths  sinee  all  equations  of  motions 
were  non-dimensionalized. 


CALCULATION  OF  X,  Y  FROM  TRANSFORMATION  EQUATION 


v2  v2co«psi2 


Figure  12.  Diagram  of  Superbloek  for  ealeulating  X  and  Y  for  vessels. 

The  tow-eonneetion  is  modeled  as  a  linear  spring  and  is  represented  by  the  final 
superbloek.  The  Kaimalino  and  Sliee  position  variables  are  used  as  inputs.  Sinee  the 
EOMs  were  derived  based  on  xg  and  yo,  the  position  variables  need  to  be  refereneed  to 
the  towing  point  eonneetion.  Equations  (34a)  through  (34d)  are  modeled  here  using 
nondimensional  length  as  the  gains.  The  spring  eonstants  Kx  and  Ky  are  eontrollable 
variables  that  are  initialized  in  the  Matlab  m-file.  These  eonstants  are  shown  in  Figure  13 
as  gains  to  the  differenee  in  Slice  /  Kaimalino  positions. 
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TRANSFER  POSITION  FROM  REFERNCE  POINT  TO  POINT  OF  TOW  CONNECTION 
AND  CALCULATION  OF  TENSION  IN  TOW  CONNECTOR 


Figure  13.  Diagram  of  Superblock  that  reference  X  and  Y  to  the  tow  connection  point 

and  calculates  tension. 
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III.  DETERMINATION  OF  HYDRODYNCAMIC  COEFFICIENTS 

A,  METHODS  FOR  DETERMINING  HYDRODYNAMIC  COEFFICIENTS 

Theoretical  predictions  of  hydrodynamic  coefficients  are  not  exact  and  still 
include  many  assumptions.  This  is  due  heavily  to  the  inability  of  theoretical  methods  to 
handle  viscous  flow  situations  that  occur  in  the  nonlinear  range.  However  correlation 
with  more  exact  experimental  methods  has  been  achieved  in  the  linear  prediction  range. 
Because  of  this  problem,  tests  have  been  done  in  order  to  provide  a  systematic  data  base 
for  predicting  hydrodynamic  coefficients.  In  Reference  [7],  the  data  obtained  during  tests 
was  coupled  with  theoretical  methods  in  order  to  guide  semi-empirical  methods  for 
calculating  hydrodynamic  coefficients.  These  tests  were  done  on  certain  geometrical 
shapes,  including  cylindrical  bodies  with  hemispherical  noses  and  cone  shaped  tails.  This 
paper  uses  the  semi-empirical  methods  derived  in  References  [6]  and  [7]. 

B,  MODEL  FOR  KAIMALINO  PODS  AND  STRUTS 

The  Kaimalino  pods  are  modeled  as  a  78  ft  long  body  of  revolution  with  a  7.4  ft 
diameter.  The  nose  of  the  pod  is  elliptic  with  a  nose  length  (Ln)  of  4.2  ft.  The  conical 
tail  of  the  pod  is  approximately  16  ft.  Attached  to  the  pod  are  the  struts  that  are  modeled 
as  24  X  5.5  ft  flat  plates.  These  struts  connect  the  pods  to  the  body  of  the  Kaimalino.  The 
profile  is  shown  below  in  Figure  14. 
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Figure  14.  Configuration  of  one  pod  with  attaehed  struts. 


The  hydrodynamie  eoeffieients  of  the  Kaimalino  pods  are  ealeulated  using  the 
DATCOM  method  for  each  pod  without  the  fins.  The  midlength  of  the  pod  is  used  as  a 
reference  point  (xm).  For  this  model,  the  reference  point  is  thirty-nine  feet  (Lb/2).  The 
coefficients  will  be  non-dimensionalized  with  respect  to  the  pod  length  (Lb).  When 
summarized  later  on,  the  coefficients  will  be  re  non-dimensionalized  with  respect  to  the 
length  of  the  Slice  model  for  consistency  in  the  towing  simulation. 

The  struts  are  modeled  as  flat  plates  fixed  fins.  The  struts  are  referenced  to  the 
center  of  the  pod  by  a  distance  Xf  that  is  the  distance  from  the  quarter  chord  point  of  the 
fin  to  the  center  of  the  pod.  After  the  individual  hydrodynamic  coefficients  for  the  struts 
have  been  obtained,  they  will  be  translated  and  non-dimensionalized  with  respect  to  the 
length  of  the  pod  and  later  re  non-dimensionalized  with  respect  to  the  SLICE  length. 

C.  HYDRODYNAMIC  COEFFICIENTS  OF  KAIMALINO  PODS 

In  order  to  maintain  consistency,  the  same  method  that  was  used  for  calculating 

the  SLICE  model  hydrodynamic  coefficients  is  used  here  in  the  calculations  for  the 

Kaimalino  model  hydrodynamic  coefficients.  The  models  have  similar  geometry.  The 
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DATCOM  method,  which  is  a  semi-emperical  method  was  used  in  the  calculation  of  all 
the  coefficients  except  for  the  acceleration  hydrodynamic  coefficients  Y'v  and  N'r  .  For 


these  two,  the  works  of  Humphrey  and  Watkinsn  [6]  was  used. 

calculated  from  the  following  equations  using  DATCOM  method 

The  coefficients  were 

o 

(Tv)pod=  '[Cla+Cno] 

1 

(35) 

,  -  ShCma 

{N'v)pod  = 

(36) 

{Tr)pod  =  ^CLl\-  — 

/  \  1  ) 

(37) 

{N'r)  pod  =  ^Cmq 

(38) 

where  Sb  is  the  maximum  cross-sectional  area  of  the  pod  and  1  is  the  total  length 
of  the  pod.  Cia  is  the  slope  of  the  lift  curve.  Cdo  is  the  drag  coefficient  of  an  elliptical 
form  at  a  zero  angle  of  attack.  Cma  is  the  slope  of  the  yawing  moment  curve.  Cmq  is  the 
slope  of  the  yawing  moment  curve  for  the  plansform  based  on  effective  area.  The  slopes 
and  values  where  calculated  with  the  following  equations 

CLa  =  2{k2-kl)—  (39) 

Sb 

where  k2  and  ki  are  the  Lamb’s  coefficients  of  inertia  for  a  prolate  ellipsoid  in 
axial  and  cross  sectional  flow  respectively.  Sv  is  the  cross  sectional  area  at  a  point  k.  Ly 
is  the  point  on  the  pod  where  flow  becomes  viscous.  The  empirical  relation  for 
calculating  k  is 

/v  =  0.378k  +  0.527k.  (40) 

Lms  is  the  distance  from  the  nose  of  the  pod  to  the  point  of  maximum  slope.  For 
this  model,  Lms  is  equal  to  Lb,  the  total  length  of  the  pod. 
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Cua  =  2 


(41) 


ki  -  k\\r  dS 

-  I  - \Xm  —  X)ClX 

y  Sbl  ){dx 


Cmq  —  Cma\ 


1- 


Xm 


f_F 

^S,b*l 


Ic  Xrn^ 

J~T 


yv 


1- 


Xm 

T 


r 


V 


Stb  *  / 


(42) 


Here  V  is  the  volume  of  the  submerged  body,  h  is  the  distanee  from  the  nose  to 
the  eenter  of  buoyaney,  and  S*  is  the  area  of  the  truneated  base  at  ly. 

The  aeeeleration  hydrodynamic  coefficients  are  calculated  as  follows  using  the 
works  of  Humphrey  and  Watkinson  [6]. 

kiV 

{Tv)pod  =  ^  (43) 

0.5/' 


{N'k)pod  =  -kbr.df  +  Y'v 

V 


Icb-Ld^ 


(44) 


where  kb  is  Lamb’s  coefficient  of  added  moment  of  inertia  for  a  prolate  ellipsoid 
body.  I’zdf  is  the  mass  moment  of  inertia  of  the  displaced  fluid  about  the  z-axis.  Table  3 
shows  the  coefficients  for  single  pods  for  the  Kaimalino  and  the  results  obtained  by 
William  Wolkerstorfer  for  a  SLICE  single  pod. 


Coefficients 

SLICE  POD 

KAIMALINO  POD 

T'v 

-0.016244 

-0.0592 

TV 

-0.002498 

-0.0014 

N'v 

-0.038754 

-0.0086 

24 


N'r 

-0.006396 

-0.0087 

T'v 

-0.044189 

-0.0407 

Tr 

0.0 

-0.0001 

N'v 

0.0 

0.0001 

N'r 

-0.001799 

-0.0045 

Table  3.  Hydrodynamic  coefficients  for  SLICE  and  Kaimalino  pods. 


D,  HYDRODYNAMIC  COEFFICIENTS  OF  KAIMALINO  STRUTS 


Theoretical  predictions  for  hydrodynamic  coefficients  for  fixed  fins  in  the  linear 
range  are  presented  in  [3].  The  process  for  theoretical  prediction  of  forces  and  moments 
in  the  nonlinear  range  is  still  an  area  under  development.  One  of  the  assumptions  that  is 
made  in  linear  theory  is  that  the  effects  of  the  ship’s  hull  on  its  fixed  fins  is  negligible. 
For  traditional  hulls,  that  assumption  is  not  really  justified  especially  in  situations  of  tight 
maneuvers.  For  the  Kaimalino  hull,  the  assumption  is  more  justified  because  of  the 
swath  hull  design.  The  hydrodynamic  coefficients  for  the  fixed  fin  are  as  follows. 


{Yv)fin  =  -A 

V 


dCi 


+  Cu 

y 


(45) 


{N\)fin  =  Y\*x'f  (46) 

{Y\)fin  =  x'f*Y\  (47) 

iN'r)fin  =  Y'v*x'/  (48) 

where  A’  is  the  fin  area  that  is  non-dimensionalized  by  the  length  times  the  draft. 
dC^ld/3  is  the  slope  of  the  lift-curve  where  |3  is  the  angle  of  attack  at  the  fin.  Cd  is  the 

drag  coefficient  at  zero  angle  of  attack.  Cd  is  small  compared  to  the  lift-curve  slope  for 
fixed  fins  therefore  it  is  ignored.  Equations  (45)  through  (52)  are  derived  from 
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summation  of  forces  and  moments  about  point  with  Xf  taken  as  the  distanee  from  the 
quarter  ehord  point  to  Xm. 


{rv)fin 


-AtA'A'j 

'ja^G  +1 


(49) 


(N'v)fin  =  Y'v*x'f 


(50) 


{Ti-)fin  =  Y'v*x'f 


(51) 


{N'r)  fin  =  Y'v*x'^  f 


(52) 


where  b  is  the  geometric  span  and  ao  is  the  geometrie  aspeet  ratio  b^/Af.  The 
effeetive  aspeet  ratio  of  a  fin  attaehed  to  ship’s  hull  beeomes  2aG. 

The  total  contribution  for  one  pod/fins  configuration  is  now  calculated  by 
simply  adding  the  contributions  from  both  forward  and  aft  fins  to  the  pod. 


(TV)  =  (T'  k) pod  +  (T'  v\fm,  aft  +  (T'  v) fin,  fwd 


(53a) 


(7V)  =  {Y'r)pod  +  {Y\)ftn,afi+{Tr)fm,fi.d 


v)  —  (  A' v) pod  +  ( A' v)  fin,  aft  +  fin,  fwd 


—  (^N'r^pod  +  (^N\\fin,aft  fin,  fwd 


(TV)  =  (T'  v) pod  +  (T'  v\fin,  aft  +  (T'  v) fin,  fwd 


(TV)  =  {T  f)pod  +  {T  r)fin,aft+{Y'  r)fin,fwd 


v)  —  (  AV)  pod  +  ( AV)  fin,  aft  +  (A'  v)  fin,  fwd 


pod  +  r\fin,  aft  +  fin,  fwd 


(53b) 

(53c) 

(53d) 

(53e) 

(530 

(53g) 

(53h) 
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E,  COMPLETE  HYDRODYNAMIC  DERIVATIVES  FOR  THE  SLICE  AND 
KAIMALINO  CONFIGURATIONS 

The  final  step  in  ealeulating  the  hydrodynamie  eoeffieients  is  to  translate  the 
eontributions  for  eaeh  pod/strut  eonfiguration.  This  is  aehieved  by  summing  forees  and 
moments  about  a  eornmon  referenee  point  on  the  vessel.  Adapting  the  equations  derived 
in  Prineiples  of  Naval  Arehiteeture  [3]  results  in  the  following  equations. 


iY\)  =  iY\)  (54) 

(TV)  =  (rv)x’/  +  (TV)  (55) 

(N\)  =  iY\)x'f  +  iN\)  (56) 

(N'r)  =  {Y\)x'f  +  (F  v)x'"  /  +  (A'v)x'/  +  {N\)  (57) 


where  x’f  is  taken  as  the  distanee  from  amidship  the  vessel  to  the  eenter  point  of 
the  pods. 

Similar  logie  is  repeated  for  the  aeeeleration  hydrodynamie  derivatives. 


(Th)  =  (Th)  (58) 

(TV)  =  (Th)x’/  +  (Th)  (59) 

{N\)  =  {Tv)x'f  +  {N'v)  (60) 

{N'  r)  =  {Y'  r)x'  f  +  iY'  v)x'"  /  +  ( A'  v)x'  f  +  iN'r)  (61) 


The  hydrodynamie  eoeffieients  were  then  non-dimensionalized  with 
respeet  to  the  eharaeteristie  length  of  the  SLICE  vessel  in  order  to  maintain  eontinuity  for 
the  towing  simulation.  Assuming  similarity  for  port  and  starboard  pod/strut 
eonfiguration,  the  eoeffieients  were  then  multiplied  by  two.  A  summary  of  both 
Kaimalino  hydrodynamie  eoeffieients  and  SLICE  eoeffieients  that  were  ealeulated  in 
previous  thesis  work  are  provided  in  the  table  below. 
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Hydrodynamic  Coefficients 

Kaimalino  (w/o  rudders) 

SLICE  (w/o  rudders) 

Tv 

-0.1183 

-0.078930 

Tr 

-0.0042 

-0.004044 

N'v 

-0.0187 

-0.016428 

N'r 

-0.0176 

-0.010332 

Tv 

-0.0184 

-0.051328 

Tv 

-0.0011 

0.005617 

N'v 

-0.0008489 

-0.001945 

N'r 

-0.0090 

-0.00564 

Table  4.  Complete  hydrodynamic  coefficients  for  SLICE  and  Kaimalino  vessels. 
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IV.  RESULTS  AND  DISCUSSION 


A,  INITIAL  MODEL  TEST 

Initially,  we  conducted  the  1®*  set  of  runs  in  order  to  validate  the  simulation 
model.  For  this  run,  the  towing  vessel  was  on  the  same  path  as  the  lead  vessel.  The 
initial  conditions  were  as  follows. 


Slice  initial  speed  =15  knots 
Kaimalino  initial  speed  =13.5  knots 
Slice  initial  yaw  angle  =  0 
Kaimalino  initial  yaw  angle  =  0  degrees 
Kaimalino  offset  from  path  =  0 


Kx  is  varied  while  Ky  is  held  constant,  and  the  longitudinal  offset  is  observed. 
Figure  15  shows  the  results  of  these  simulation  runs.  The  trend  shows  a  decreasing  offset 
for  higher  values  of  stiffness  in  the  Kx  parameter  as  expected.  The  non-dimensional 
value  Kx  value  of  0.06  corresponds  to  a  near  rigid  towing  connection  where  there  is  very 
little  deformation  based  on  the  stiffness. 
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LONGITUDINAL  OFFSET 


Figure  15.  Longitudinal  Offset  observation  while  varying  Kx  for  3rd  set  of  simulation 

runs  initial  eonditions. 

B.  DIRECTIONAL  STABILITY 

Various  simulation  runs  were  conducted  in  order  to  analyze  the  directional 
stability  of  the  towing  model.  There  are  different  kinds  of  motion  stability  when  dealing 
with  surface  vessels.  In  general,  motion  stability  is  the  ability  of  the  vessel  to  retain  as 
much  of  the  initial  state  of  equilibrium  in  its  final  path  after  some  disturbance  has 
occurred.  For  directional  stability,  the  vessel’s  final  path  will  retain  the  straight  line 
motion  it  initially  had  as  well  as  its  initial  direction.  A  vessel  may  oscillate  after  the 
disturbance  or  may  transition  smoothly  back  to  its  original  state.  The  latter  is  obviously 
preferred  since  it  corresponds  to  a  more  stable  vessel. 

In  the  simulation  runs  conducted  for  this  paper,  two  sets  of  initial  conditions  were 
applied.  For  these  runs,  the  vessel  being  towed  has  experienced  a  disturbance  and  is 
offset  from  the  initial  path  set  by  the  lead  vessel.  The  lateral  offset  is  observed  to  study  if 
the  vessel  is  directionally  stable  or  not.  For  each  of  these  sets,  different  combinations  of 
stiffness  coefficients  are  applied.  For  the  L'  set  of  runs,  the  initial  conditions  are  as 
follows. 


Slice  initial  speed  =15  knots 
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Kaimalino  initial  speed  =13.5  knots 


Slice  initial  yaw  angle  =  0 


Kaimalino  initial  yaw  angle  =  0 


Kaimalino  Offset  from  path  =  0.1  (non-dimens ionalized  to  SLICE  length) 


Figures  16  and  17  show  the  results  of  varying  Ky  while  holding  Kx  constant. 


LATERAL  OFFSET 


Figure  16.  Fateral  offset  for  varying  Ky  for  F*  set  initial  conditions. 


KAIMALINO  YAW  ANGLE 


Figure  17.  Kaimalino  yaw  angle  for  varying  Ky  for  F*  set  initial  conditions. 
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From  the  figures  it  is  evident  that  as  Ky  increases,  the  lateral  offset  and  the  yaw 
angle  oscillate  more.  This  indicates  that  the  system  will  eventually  go  unstable  if  Ky  is 
increased  too  much  while  holding  Kx  constant.  The  oscillation  is  around  the  zero  offset 
which  is  the  critical  point  to  be  used  in  an  eigenvalue  analysis. 

The  next  two  figures  show  the  results  of  varying  Kx  while  holding  Ky  constant. 


LATERAL  OFFSET 


Figure  18.  Lateral  offset  for  varying  Kx  for  set  initial  conditions. 


KAIMALINO  YAW  ANGLE 


Figure  19.  Kaimalino  yaw  angle  for  varying  Kx  for  L*  set  initial  conditions. 
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Now  for  an  increasing  Kx  value,  there  is  significantly  less  oscillation  and  the 
systems  is  more  stable  at  the  higher  values. 

The  last  results  for  the  set  of  initial  conditions  are  for  varying  Kx  and  Ky  at  the 
same  time. 


LATERAL  OFFSET 


Figure  20.  Lateral  offset  for  varying  Kx  and  Ky  for  L'  set  initial  conditions. 


KAIMALINO  YAW  ANGLE 


Figure  21 .  Kaimalino  yaw  angle  for  varying  Kx  and  Ky  for  set  initial  eonditions. 
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The  figures  show  that  when  inereasing  both  Kx  and  Ky,  the  direetional  stability 
deereases.  However,  inereasing  Kx,  slows  down  the  oscillations  that  were  occurring  in 
the  case  of  increasing  Ky  only. 

For  the  2nd  set  of  runs,  the  initial  conditions  are  as  follows. 


Slice  initial  speed  =15  knots 
Kaimalino  initial  speed  =13.5  knots 
Slice  initial  yaw  angle  =  0 
Kaimalino  initial  yaw  angle  =10  degrees 
Kaimalino  Offset  from  path  =  0 


Figures  22  through  27  show  the  results  of  different  combinations  of  stiffness 
coefficients  for  the  2"‘^  set  of  initial  conditions.  The  trend  for  the  results  is  similar  to  what 
was  observed  in  the  set  of  simulation  runs.  For  the  Kaimalino  yaw  angles,  the  results 
show  that  the  angles  return  to  the  critical  point  (0  offset)  either  smoothly  or  with 
oscillation.  Either  case  the  system  has  some  measure  of  stability  for  the  different 
stiffness  coefficients. 


LATERAL  OFFSET 
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OFFSET 


Figure  22.  Lateral  offset  for  varying  Ky  for  2nd  set  initial  conditions. 


KAIMALINO  YAW  ANGLE 


Figure  23.  Kaimalino  yaw  angle  for  varying  Ky  for  2nd  set  initial  conditions. 


LATERAL  OFFSET 


Figure  24.  Lateral  offset  for  varying  Kx  for  2nd  set  initial  conditions. 
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OFFSET 


KAJMALINO  YAW  ANGLE 


Figure  25.  Kaimalino  yaw  angle  for  varying  Kx  for  2nd  set  initial  eonditions. 


LATERAL  OFFSET 


Figure  26.  Lateral  offset  for  varying  Kx  and  Ky  for  2nd  set  initial  conditions. 
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KAJMALINO  YAW  ANGLE 


Figure  27.  Kaimalino  yaw  angle  for  varying  Kx  and  Ky  for  2nd  set  initial  conditions. 
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V.  CONCLUSIONS 


A  simulation  model  has  been  developed  that  represents  the  eoupled  dynamies  of 
two  vessels  in  a  towing  operation.  This  nonlinear  time-dependent  simulation  model  is 
useful  for  eourse  stability  analysis  of  more  sophistieated  towing  systems.  The  model 
takes  into  aeeount  the  effeets  of  towline  dynamies  through  a  set  of  generalized  stiffness 
eoeffleients.  Based  on  the  results  obtained  from  the  simulations  runs  presented  in  this 
paper,  it  is  feasible  to  have  a  near  rigid  towing  eonneetion  for  the  Sliee  towing  the 
Kaimalino  in  elose  proximity.  However,  the  results  seem  to  indieate  that  the  need  for 
some  flexibility  in  the  longitudinal  direetion  for  the  stiffness  or  even  a  eontrolled 
attaehment  point  whieh  would  allow  for  longitudinal  motion  of  the  towing  eonneetor. 

As  a  reeommendation  for  further  researeh,  a  study  of  the  eigenvalues  using  a 
perturbation  analysis  near  the  eritieal  points  should  be  eondueted  in  order  to  find  the 
range  of  Kx  and  Ky  values  that  make  the  system  stable  or  unstable 
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APPENDIX  A.  MATLAB  PROGRAM  FOR  CALCULATION  OF 
HYDRODYNAMIC  COEFFICIENTS  FOR  SINGLE  POD/STRUT 
CONFIGURATION  AND  COMPLETE  HYDRODYNAMIC 
COEFFICIENTS. 


%LT  Garrett  Jones 
%Thesis  Work 

%Kaimalino  Hydrodynamie  Coeffieients 

%Following  programs  are  developed  to  ealeulation  Hydrodynamie  Coeffieients  (HCs) 
%for  Kaimalino  using  the  DATCOM  method  and  works  of  Humphrey  and  Watkinson. 
%This  method  was  used  by  William  Wolkerstorger  in  ealeulating  HCs  for  the  SLICE 

%Single  Pod  ealeulations 

%Initial  Values 
gamma=62.4;  %lbf/fG3 
g=32.2;  %ft/see^2 
rho  =gamma/g; 
r=4;  %  radius  of  pod 
lb=78;  %  base  length 

lms=78;  %  length  to  point  of  maximum  slope 

ln=4;  %  nose  length 

lf=62;  %  length  of  forebody 

lo=lb-lf;  %  length  of  oonieal  tail 

lpmb=lf-ln;  %  length  of  parallel  midbody 

xm=39;  %  referenee  point  for  ealeulations 

d=7.4; 

%  ealeulations  of  seetional  areas,  eurve  slopes,  and  eoeffieients  needed 
%  needed  for  determining  HCs 

lv=0.378*lb+0.527*lms;  %  length  where  flow  beeomes  viseous 
Sv=pi*(r*(lb-lv)/lo)^2; 
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Sb=pi*r^2; 

k2=.963269; 

kl=.019066; 

kb  =0.892252; 

Cla=2*(k2-kl)*(Sv/Sb); 

Cd=0.29; 

INTl=(-0.5*pi*d^2/ln^2)*((l/3)*ln^3-(l/2)*ln^3+xm*ln^2-(l/2)*xm*ln^2); 

INT2=(.5*pi*d^2/lc^2)*[(-.5*lv^2*lc+(l/3)*lv^3-.5*lv^2*lf+xm*lc*lv-.5*xm*lv^2+... 

xm*lf*lv)-(-.5*lf^2*lc+(l/3)*lf^3-.5*lf^3+xm*lc*lf-.5*xm*lf^2+xm*lf^2)]; 

Cma=(2*(k2-kl)/(Sb*lb))*(INTl+INT2); 

Stb=(d/lc)*(lc*lb-(l/2)*lb^2+lPlb-lc*lv+(l/2)*lv^2-lPlv); 

V=(pi*d^2/4)*(ln-ln/3+lpmb+lb^3/(3*lc^2)-lb^2*lf/lc^2-lP3/(3*lc^2)+... 

lb*lf^2/lc^2); 

lcb=[(pi*d^2/4)*((2/3)*ln^2-ln^2/4+0.5*lpmb^2+lbM/(12*lc^2)-lb^2*... 

lf^2/(2*lc^2)-lfM/(4*lc^2)+2*lb*lP3/(3*lc^2))]A^; 

leg  =lcb; 

Cmq=Cma*  [((1  -xm/lb)^2-(  V/(Stb  *  lb))  *  ((lcb/lb)-(xm/lb)))/ (( 1  -xm/lb)-(  V/. . . 

(Stb=^lb)))]; 

C  lq=Cla  *  ( 1  -xm/lb) ; 

lydf  =(rho*pi*d^2/4)*[(3*ln^3/10)-(5*xm*ln^2/6)+(2*xm^2*ln/3)+xm^2*lpmb-... 
xm*lpmb^2+(lpmb^3/3)+(l/lc^2)*((lb^5/30-xm*lbM/6+xm^2*lb^3/3)-... 
(lf^5/5-xm*lfM/2-lb*lfM/2+xm^2*lf^3/3+lb^2*ir3/3+4*xm*lb*lf^3/3-... 
xm^2*lb*lf^2-xm*lb^2*lf^2+xm^2*lb^2*lf))]; 


%  Calculations  of  HCs 
Yv=-(Sb/lb^2)*(Cla+Cd); 

Yr  =  -(Sb*Clq)/lb^2; 

Nv=(Sb/lb^2)*Cma; 

Nr  =  -(Sb*Cmq)/lb^2; 

Yvdot  =  -(k2*V)/(0.5*lb^3); 

Nrdot  =  (-kb*Iydf/(0.5*rho*lb^5))+Yvdot*((lcb-lcg)/lb); 
HCP(l)=Yv; 

HCP(2)=Yr; 
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HCP(3)=Nv; 

HCP(4)=Nr; 

HCP(5)=Yvdot; 

HCP(6)=0; 

HCP(7)=0; 

HCP(8)=Nrdot; 
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%LT  Garrett  Jones 
%Thesis  Work 

%Kaimalino  Hydrodynamie  Coeffieients 

%  Struts  Caleulations 

% - 

%  Initial  Values 

LB  =  78;  %length  of  pod 

Lbp  =  81.8;  %length  between  perpendieulars 

T=14;  %  draft 

xm=39; 

xf=28.8;  %distanee  from  xm  to  quarter  chord  fwd  fin 
xf2=-29.2;  %distance  from  xm  to  quarter  chord  aft  fin 
b=5.5; 
c=24; 


%  Calculation  of  areas,  aspect  ratios,  and  curve  slopes  needed  for  determing  HCs 

A=b*c; 

ag=2*(b/c); 

a=ag; 

Cla=(0.9*pi*2*a)/(sqrt(a^2+4)+1.8); 

Cmcl=  1  /2-((  1 . 1 1  *  sqrt(ag^2+4)+2)/ (4  *ag+2)); 

Cmc4=(.25-Cmcl)*Cla; 

Xf=xf; 

Xf2=xf2; 

%  HCs  calculations  for  single  strut 
Yv=-abs(A*Cla/(LB*T)); 

Yr=(Xf/LB)*Yv; 

Nv=Yr; 

Nr=(Xf/LB)^2*Yv; 

Y  vdot=-(2  *pi*  c  *b^2)/ ((LB^2  *  T)  *sqrt(ag^2+ 1 )); 
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Y  rdot=-(Xf/LB)*  Yvdot; 

Nvdot=(Xf/LB)*  Y  vdot; 

Nrdot=(Xf/LB)^2*Yvdot; 

Yv2=-abs(A*Cla/(LB*T)); 

Yr2=(Xf2/LB)*Yv2; 

Nv2=Yr2; 

Nr2=(Xf2/LB)^2*Yv2; 

Yvdot2=-(2*pi*c*b^2)/((LB^2*T)*  sqrt(ag^2+ 1 )) ; 
Yrdot2=-(Xf2/LB)*Yvdot2; 
Nvdot2=(Xf2/LB)*Yvdot; 
Nrdot2=(Xf2/LB)^2*Yvdot2; 


HCS(1)  =  Yv; 
HCS(2)  =  Yr; 
HCS(3)  =  Nv; 
HCS(4)  =  Nr; 
HCS(5)  =  Yvdot; 
HCS(6)  =  Yrdot; 
HCS(7)  =  Nvdot; 
HCS(8)  =  Nrdot; 


HCS2(1)  =  Yv2; 
HCS2(2)  =  Yr2; 
HCS2(3)  =  Nv2; 
HCS2(4)=Nr2; 
HCS2(5)  =  Yvdot2; 
HCS2(6)  =  Yrdot2; 
HCS2(7)=Nvdot2; 
HCS2(8)  =  Nrdot2; 


% - Non-dimensionalize  and  combining  strut  and  pod  coefficients - 

HC(1)=  HCP(1)(LB/Lbp)^2  +  2HCS(1)*(LB*T)/Lbp^2+  2HCS2(1)(LB*T)/Lbp^2; 
HC(2)=  HCP(2)(LB/Lbp)^3  +  2HCS(2)(LB^2*T)/Lbp^3+  2HCS2(2)(LB^2*T)/Lbp^3; 
HC(3)=  HCP(3)(LB/Lbp)^3  +  2HCS(3)(LB^2*T)/Lbp^3+  2HCS2(3)(LB^2*T)/Lbp^3; 
HC(4)=  HCP(4)(LB/Lbp)M  +  2HCS(4)(LB^3*T)/LbpM+  2HCS2(4)(LB^3*T)/Lbp^4; 
HC(5)=  HCP(5)(LB/Lbp)^3  +  2HCS(5)(LB^2*T)/Lbp^3+  2HCS2(5)(LB^2*T)/Lbp^3; 
HC(6)=  HCP(6)(LB/Lbp)M  +  2HCS(6)(LB^3*T)/LbpM+  2HCS2(6)(LB^3*T)/Lbp^4; 
HC(7)=  HCP(7)(LB/Lbp)M  +  2HCS(7)(LB^3*T)/LbpM+  2HCS2(7)(LB^3*T)/Lbp^4; 
HC(8)=  HCP(8)(LB/Lbp)^5  +  2HCS(8)(LBM*T)/Lbp^5+  2HCS2(8)(LB^4*T)/Lbp^5; 


%■ 


■Port  and  Stbd  HCs 
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Yvhc  =  2*HC(l) 
Yrhc  =  2*HC(2) 
Nvhc  =  2*HC(3) 
Nrhc  =  2*HC(4) 
Yvdothc  =  2*HC(5) 
Yrdothc  =  2*HC(6) 
Nvdothc  =  2*HC(7) 
Nrdothc  =  2*HC(8) 
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